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Abstract 

We consider minimal modifications to tribimaximal (TBM) mixing matrix which accommodate 
non-zero mixing angle #13 an d CP violation. We derive four possible forms for the minimal modifi¬ 
cations to TBM mixing in a model with A 4 flavor symmetry by incorporating symmetry breaking 
terms appropriately. We show how possible values of the Dirac-type CP phase dp can be predicted 
with regards to two neutrino mixing angles in the standard parametrization of the neutrino mixing 
matrix. Carrying out numerical analysis based on the recent updated experimental results for neu¬ 
trino mixing angles, we predict the values of the CP phase for all possible cases. We also confront 
our predictions of the CP phase with the updated fit. 
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I. INTRODUCTION 


Establishing leptonic CP violation (LCPV) is one of the most challenging tasks in future 
neutrino experiments [I]. The relatively large value of the reactor mixing angle measured 
with a high precision in neutrino epxeriments [2] has opened up a wide range of possibilities 
to explore CP violation in the lepton sector. The LCPV can be induced by the PMNS 
neutrino mixing matrix [3] which contains, in addition to the three angles, a Dirac type CP 
violating phase in general as it exists in the quark sector, and two extra phases if neutrinos 
are Majorana particles. Although we do not yet have compelling evidence for LCPV, the 
current global fit to available neutrino data indicates nontrivial values of the Dirac-type CP 
phase HIE]. In this situation, it must deserve to predict possible size of LCPV detectable 
through neutrino oscillations. From the point of view of calculability, much attention has 
been paid to the prediction of the Dirac type LCPV phase with regards to some observables 
[Sj. Recently, it has been shown [7] that Dirac-type leptonic CP phase can be particularly 
predictable in terms of neutrino mixing angles in the standard parameterization of PMNS 
mixing matrix [8]. 

Before the measurements of the reactor mixing angle, the fit to neutrino data was con¬ 
sistent with the so-called tribimaximal (TBM) neutrino mixing matrix, which is 

theoretically well motivated flavor mixing pattern pjj. However, it should be modified to 
accommodate non-zero reactor mixing angle as well as CP violation. Although the current 
neutrino data rule out the exact TBM mixing pattern, it can be regarded as leading or¬ 
der approximation. Among various possible modification to Ltbm, as discussed in [7], the 
minimal modihcaton is useful to predict Dirac type CP phase. The minimal modification is 
to multiply t/™ M by a rotation matrix in the (i , j ) plane with an angle 6 and a CP phase 
£, Uij(6, £), whose form is given either U}j(Q,£)Uj BU or U^ nM U l: j(9, £) [ID] . Among them, 
Ul 3 (9, £,)Uq BM and Uq BM Ui 2 (0, £) are ruled out because they lead to zero reactor mixing 
angle. So, all possible forms of minimal modification to TBM mixing matrix are as follows: 


UTBUi 


V = 


o U 23 (9,0 (Case-A), 

U™U 13 (6,0 (Case-B), 

A(«.e)C„ TBM (Case-C), 
l c4(»,f)C/ 0 ™ (Case-D). 

While the study in |7j has not accounted for the origin of such modification to U ( 


( 1 ) 


TBM 
0 ) 
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in this paper, we first study how such a minimally modified TBM mixing pattern can 
be achieved in a neutrino model with A 4 flavor symmetry by incorporating A 4 symmetry 
breaking terms appropriately. Then, following [7], we investigate how the Dirac type CP 
phase can be predicted based on the updated fit results for neutrino mixing angles [5]- As 
shown later, comparing with the results obtained in [Tj, the Dirac type CP phase predicted 
based on the updated fit results has different implication particularly at la C.L. 

II. MINIMAL MODIFICATIONS TO TRI-BIMAXIMAL MIXING IN A 4 SYM¬ 
METRIC MODEL 

In mi, an A 4 symmetric model for neutrino masses and mixing has been proposed to 
accommodate non-zero mixing angle 0 13 on top of TBM mixing. Based on the A 4 symmetric 
model, we study how the forms given in Eq.Q can be derived by incorporating appropriate 
A 4 symmetry breaking terms. 


Case-A 


As proposed in |llj . A 4 flavor symmetry allows the charged-lepton mass matrix to be 
diagonalized by the Cabibbo-Wolfenstein matrix ra, 


1 1 1 1 


Ucw — 


a/3 


1 U! U ) 2 
1 Ul 2 00 


( 2 ) 


where u = e 2m / 3 , with three independent eigenvalues, rrie^m^^mr. This can be realized by 
the lepton assignments, U = M) ~ 3, i\ ~ 1 , 1% ~ v, 11 ~ i" with 3 Higgs doublets 
= {<,~ 3. Introducing 6 heavy A 4 Higgs singlets and triplet: 


m ~ i, V2 ~ I', vs ~ I", Vi(= 4,5,6) ~ 3 , 

where r/, = (7 7 7 ^, 77 A, 77? ), one can obtain the neutrino mass matrix in the A 4 basis [T]L| 


(3) 


/ 


M v — 


a + b + c f e 

f a + cob + co 2 c d 

e d a + uj 2 b + 00 c 


\ 


(4) 
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where a comes from < rf( >, b from < 77 ° >, c from < 77 ° >, d from < 77 ° >, e from < 77 ° >, 
/ from < 7 ]^ >. To achieve TBM mixing pattern of the neutrino mixing matrix, A 4 flavor 
symmetry should be broken to Z 2 in such a way that b = c and e = f = 0. Then, the 
neutrino mass matrix in the flavor basis where the charged lepton mass matrix is diagonal 
is given by 


( a + (2d/3) b - (d/3) b - (d/3) ^ 


/or( e )M> T ) — rrt jur tj* _ 
lvl v — u CW nl v U CW — 


b - (d/3) b + (2d/3) a - (d/3) , (5) 

^ 6 — (d/3) a— (d/3) 6+(2d/3)y 

which is diagonalized by the TBM mixing matrix C/’o' BM . To achieve non-zero mixing angle 
$13 so as to accommodate neutrino data from reactor experiments, we take b = c and 
e = —/ = e 7 ^ 0 in Eq.Q, and then the neutrino mass matrix in the A 4 basis is given by 


^ a + 2 b e -e ^ 


M„ = 


\ 


a — b d 
d a — b 


( 6 ) 


In the flavor basis, the neutrino mass matrix can be rewritten as 


M (e,M ’ T) = 


/ 


+ 


v/3 


1 a + (2d/3) b - (d/3) b - (d/3) ^ 

6 —(d/3) 6+(2d/3) a— (d/3) 

^ 6 —(d/3) a — (d/3) 6+(2d/3)y 

Rotating the mass matrix given in Eq. ([7j) by TBM mixing matrix, we get 

/ ^ u 1 a n n \ 


V 


0 -e e 
-e — 2 e 0 
e 0 2 e 


(7) 


a — b + d 0 0 

0 a + 2b X 
0 X b — a + d 


( 8 ) 


where X = \[2ie and non-zero entries are complex in general. It can be easily shown that 
the mass matrix given by Eq. (J7]) can be diagonalized by 


/ 


V' = U™ 


1 0 0 
0 cos 6 — sin 

0 sin 9e 1 ^ cos 6 


\ 


■Ph 


(9) 


where Pp = Diag[e l/31 , e 1 ^ 2 , e^ 3 ]. 
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Now, let us check testability of the cases in this neutrino model by taking into account the 
sum-rules among the light neutrino masses da. In the leading order, the mass eigenvalues 
are given by rn® = a — b + d, m® = a + 2 b, = b — a + d, and thus we get the mass sum 

rules 


m® = m 2 + 771 ij f° r a = 0, 

m® = 2 m® + m 3 i f° r b = 0. (10) 

Inclusing the perturbation given by the second matrix in Eq. (Jr]), we get the following sum 
rule, 


ifi 2 + m 3 = hi ° + m®, 


( 11 ) 


where m 2 = m 2 ^ m 3 = 7 / 736 ^, rh® = m®e rh® = m®e^ with m^= 1 , 2 , 3 ) representing the 


mass eigenvalues obtained by diagonalizing the mass matrix Eq.(J 8 j). Plugging Eq.(10) into 
Eq.gg, we can get the following sum rules for £ = 0 , n, 2 r r, 


m\ + m 2 — 7773 = 2 $ 777 - 2 ) for a = 0 , 

2 7772 + 7773 — 777! = $7772, f° r $ = 0 , ( 12 ) 


where $? 7 i 2 = m 2 — m® and we have used mi = m®. The sum rules for £ = 7 t/ 2 , ( 37 t/ 2 ) are 


777! + 777 2 = m 3 , for O = 0, 
2??72 + m 3 — mi = 3 $ 7772 , for b = 0. 


( 13 ) 


B. Case-B 


To realize the case B, we add the breaking terms 5M V to M u in the A 4 basis, which is 
given by 


where g 


5M V 


f g + h 0 0 ^ 

0 ujg + uj 2 h 0 
^ 0 0 c o 2 g + uh y 


—h and A = \fiig. 


^0 0 0 ^ 
0 A 0 

^0 0 -A y 


( 14 ) 
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Then, the matrix given in Eq. (14) becomes in the flavor basis as follows: 


/ 


\ 


o g -g 
9-9 0 
\-g 0 g 

Then, the mass matrix M v + 5M U can be diagonalized by 


( 15 ) 


V = U, 


TBM 

0 


/ 


\ 


cos 9 0 — sin 9e ^ 
0 1 0 

sin 9e l( - 0 cos 9 


■Pr. 


(16) 


For the case B, the sum rules at the leading order are the same as Eq.(10). Including the 


perturbation given by the second matrix in Eq. (15), we get the following sum rule 


nil + fh-3 = m 5 + mi], 


(17) 


where rhi = m\e m 3 = in® = m\e m\ = m^e^, with mj( = i ) 2 , 3 ) representing the 


mass eigenvalues obtained by diagonalizing the mass matrix M v + 5M V . Plugging Eq.(10) 
into Eq. ([17]) , we can get the following sum rules for £ = 0, it, 27 r, 


m 3 — 777,2 — mi = 2 5777 , 3 , for a = 0, 
m 3 + 2 m 2 — mi = 2 67773 , for 6 = 0, 


(18) 


where 67773 = m 3 — 7773 and we have used m 2 = m®. The sum rules for £ = 7 t/ 2 , ( 37 t/ 2 ) are 


777 1 — 7773 + 7772 = 0, for 0 = 0, 
mi — 2 m 2 — m 3 = 0 , for 6 = 0 . 


(19) 


C. Case-C 


The case C can be realized by adding the A 4 breaking term SM t to the charged lepton 
mass matrix M; : 


SM t 


* giv 1 g 2 v 1 0 ^ 
£1^2 .92^2 0 
93^3 0 j 


( 20 ) 
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Taking v\ = v 2 = v 3 and g\ — g 2 — g, the matrix given by (20) becomes 


SM t = 


( gv gv 0 ^ 
guv gv 0 
y gu 2 v gv 0 J 




— Ucw 


g o o 
ooo 




V. 


( 21 ) 


Due to the addition of 8 M t , the PMNS mixing matrix should be changed to Ul 2 (9, C,)Uq BM P p. 


For Case C, the sum rules are given by Eq.(10). 


D. Case-D 


Similarity, the case D can be achieved by adding the following matrix 5Mi to the charged 
lepton mass matrix Mp 


8 M l = 


^ givi 0 g 2 v 1 ^ 
g\u 2 v 2 0 g 2 v 2 
y gxuv 3 0 g 3 v 3 J 


( 22 ) 


Taking V\ — v 2 — v 3 and ,9i = g 2 — g, the matrix given in (22) becomes 


5M, = 


^ gv 0 gv ^ 
gu 2 v 0 gv 
^ guv 0 gvJ 


= U, 


cw 


^0 0 g^ 
0 0 0 
\g o o j 


V3■ 


v. 


(23) 


The addition of 5Mi causes the PMNS mixing matrix changed to U{ 3 (9, £)Uq M Pp 


For Case D, the sum rules are given by Eq.(10). 


III. PREDICTIONS OF DIRAC-TYPE CP PHASE 

Now, let us review how to predict Dirac-type CP phase in PMNS mixing matrix with 
regards to neutrino mixing angles presented in (7]J. Multiplying V given in Eq.Q by phase 
matrices P Q and Pp that can be arisen from the charged lepton sector and neutrino sector, 
respectively, we can equate it with the standard parameterization of the PMN mixing matrix 
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as follows: 


P„-V-P a = U ST = u, 


rPMNS 


■p* 


( 


C12C13 


S12C13 


3 —iSo 


\ 


-S12C23 - Ci 2 S23Sl 3 e*' 5D C12C23 - Sl2S23Sl3e*' 5D 


p*. 


S I 3 e 
S23C13 

S12S23 — C 12 C 2 3S 13 e'''' i ' — C12S23 — Si 2 C23Si3e"“ u C23C13 

The equivalence between both parameterizations dictates the following relations, 


V 


0 i5n 


0 i&D 


(24) 


/ 


Vij e i ^ i+ ^ = Uf 7 = . 


v 




)ij t 


(25) 


A. Case A and B 


Applying |Vi 3 | = |and IV 11 /V 12 I = C// /f/j 3 / |, we obtain the relations 


ST /rrSTi 


sin 2 6 = 3s 2 3 , 


cos 2 6 = 


(26) 


S 12 — 


2 tan 2 032 (CaseA), 

| cot 2 0i2 (CaseB), 

which lead to the relation between the solar and reactor mixing angles, 

1 - 3 ( 14 ^) ( CaSeA )> 

a(I=i£) (CaseB). 

Those relations indicate that non-zero values of sf 3 lead to s 2 2 < 1/3 for Case A and 
s 2 2 > 1/3 for Case B . From IV^/V^I = |tl| 3 T /t/| 3 r |, we also get the relations 

(CaseA), 


(27) 


COSIjI = l 


q 2 

b 13 


(CaseB). 


(28) 


si3\/2-3 sf 

Now, we demonstrate how to derive So in terms of neutrino mixing angles in the standard 
parametrization. From the components of the neutrino mixing matrix for Case-A, we see 
that 

^23 + V33 _ 113 
^22 + ^32 1 12 


(29) 


From the relation (25), we get the relations 

UfJ _ U 23 + Ul7e~^ a 3 -“ 2 ) 
C 2 S 9 T + ’ 


Ilf? 

: U] ... 

eg T Vx 


(30) 


( 31 ) 















(32) 


Since V 2 i = VI 


'31> 


r/ ST 

,, 1 ( 0:3 — 02 ) _ ^31 

' u*f 


Plugging Eq.(32) into Eq.(30), we finally obtain the relation 


Uf 2 Ul?Ul? + Ul?Ul?' 


(33) 


Notice that the Majorana phases in Eq.(33) are cancelled. Presenting Uf^ explicitly in terms 


of the neutrino mixing angles as well as Sp, we get the equation for 5 D as 


cos Sd = 


-1 


1 - 5 * 2 3 


2 tan 2 d 2 3 si 3-\/2 — 6 s 


2 

13 


(34) 


Notice that the imaginary part in Eq. (33) is automatically cancelled. 


Similarity, we get the relation for Case B, 

c 1 

COS On = 


2 - 44 


2 tan 2 9- 


23 S13 




3 s 2 

OS 13 


(35) 


B. Case C and D 


Applying |Vi 3 | = |[/ 1 s 3 t | and \V 2i /V^\ = \U^/U^\, we obtain the relations 


sin 2 6 = 2 s 2 3 , 


cos 2 6 = 


tan 2 $23 (CaseC), 
cot 2 6*23 (CaseD), 

which lead to the relation between the atmospheric and reactor mixing angles, 


(36) 


1 - 


s 23 — 


2 (i=i£) (CaseC), 
2 aifj (CaseD). 


(37) 


Those relations indicate that non-zero values of s 2 3 lead to S 23 < 1/2 for Case C and 
s 2 3 > 1/2 for Case D . From \V U /V 12 \ = \Uf?/U??\, we also get the relation 

1 


cos ?/1 = 


3c 2 ~ 2 


13 S 12 


2313 -^/ 2 ” 


4s 13 


(38) 


We note that both cases lead to the same relation for | cos rj |. 


Following the same procedures for obtaining Eqs.(34, 35), we get the relations 


giazilzMaKlzMa) 


COS 5n = < 6si2C i 2S i3\ /l - 2 4 

' (1—SsfalC 1 —3s? 3 )—s 


(CaseC), 

(CaseD). 


(39) 

_ ’i 3 (PoDoni 

6si2Ci2Sl3\/l-2sf 

Sustituting experimental values for neutrino mixing angles into Eqs.(34j35|39), we can esti¬ 
mate the values of Sp for each cases. 
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IV. NUMERICAL RESULTS 


For our numerical analysis, we take the current experimental data for three neutrino 
mixing angles as inputs, which are given at lcr — 3cr C.L., as presented in Ref. [5J. This 
analysis is, in fact, to update the numerical results for the prediction of So given in 0 by 
taking the new fit to the data [5]. However, as shown later, the results based on lcr data 
is completely different from those in [7j. Here, we perform numerical analysis and present 
results only for normal hierarchical neutrino mass spectrum. It is straight-forward to get 
numerical results for the inverted hierarchical case, and we anticipate that the conclusion is 
not severly changed in the inverted hierarchical case. Using experimental results for three 
neutrino mixing angles, we first check if the relations Eqs.( 27p7 ) hold and then estimate 
the values of So in terms of neutrino mixing angles for those four cases. 



In 


§.21 0.28 0.29 0.30 0.31 0.32 0.33 0 . 34 ° 

,2 

S 12 


FIG. 1. Prediction of So in terms of sf 2 for Case C based on lcr experimental data. 
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A. Results for 1 a C.L. 


Plugging the experimental data for s 2 3 at lcr C.L. into Eqs.( 27p7 ), we predict the values 
of the mixing parameters sf 2 (Case A and B) and s | 3 (Case C and D) as follows: 


S 12 ~ 


s 23 — 


(40) 

(41) 


0.318 -0.319(CaseA), 

0.340 - 0.341 (CaseB), 

0.488 - 0.489(CaseC), 

0.510-0.511(CaseD). 

We see that s 2 2 and s | 3 are very narrowly determined for the la region of s 2 3 . Comparing 
the experimental values of s\ 2 and s 23 with the above predictions, we see that only Case C 
is consistent with experimental results at la C.L. 

In Fig. [lj we show the prediction of 5d in terms of s 2 2 based on the experimental data 
at la C.L. The upper curve in Fig. [l] indicates 1.327T < S D < 1.527T which is consistent with 
the result of fit for CP phase ( 1 . 37 T < 5d < 1.927t) shown in [5j. 


B. Results for 3 a C.L. 


Plugging the experimental data for sin 2 $13 at 3a C.L. into Eqs.( |27|[37| , we predict the 
values of the mixing parameters sin 2 #12 (Case A and B) and sin 2 #23 (Case C and D) as 
follows: 


S 12 — 


S 23 — 


0.316 - 0.321 (CaseA), 

0.340 -0.342(CaseB), 

(42) 

0.487- 0.491(CaseC), 

0.509 - 0.513(CaseD). 

(43) 


Comparing the experimental values of s 2 2 and s 23 with the above predictions, we see that 
they are all consistent with experimental results at 3a C.L. I 11 particular, the prediction of 


’12 


for Case B prefers to nearly upper limit of 3a allowed region. 

Fig. [ 2 ] shows the predictions of 8 d in terms of s 23 ((a): Cases A and B) and s 2 2 ((b): 
Cases C and D) based on the corresponding experimental data given at 3a C.L. Regions 
surrounded by blue and red lines correspond to Cases (A, C) and (B, D), respectively. The 
width of each bands implies the variation of the other mixing angles, s 2 2 (Cases A and B) 
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and s| 3 (Cases C and D). We see that almost maximal 8 jj ~ 7r/2,37r/2 can be achieved by 
s| 3 ~ 0.5 for Cases-A, B and by s 32 ~ 0.325 for Cases-C,D . The values around 3ir/2 is 
consistent with the current fit of the Dirac type CP phase [5]. 

Comparing those results with the corresponding ones presented in [ 7 ], we see that the 
shapes of the curves in each cases are nearly unchanged, but the widths of each bands get 
much narrower. The allowed regions of sf 2 above 0.344 for Cases C and D are excluded in 
the updated analysis. 


(a) 




<? §21 0.28 0.29 0.30 0.31 0.32 0.33 0 . 34 ° 


*12 


FIG. 2. Prediction of 5d in terms of (a) s 23 for Cases A and B, and (b) sf 2 for Cases C and D 
based on 3cr experimental data. Regions surrounded by blue (red) lines correspond to Cases A, C 
and (B, D). 


V. CONCLUSION 

As a summary, we have considered how non-zero mixing angle $i 3 and CP violation can 
be acommodated in a model with A 4 flavor symmetry by incorporating symmetry breaking 
terms appropriately. The four possible forms of neutrino mixing matrix we considered are 
minimal modifications to TBM mixing matrix and factorized by TBM mixing form and 
an unitary mixing matrix with an angle and a CP phase corresponding to a rotation in a 
plane. We have shown that possible size of the Dirac-type CP phase 5n can be predicted 
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with regards to two neutrino mixing angles in the standard parametrization of the neutrino 
mixing matrix. This has been achieved by equating one of minimally modified TBM mixing 
matrix with the standard parametrization of the PMNS one. Based on the current fit results 
for the neutrino mixing angles and CP phase, we have seen that the neutrino mixing matrix 
corresponding to Case C is consistent with the current fit data at lcr C.L. whereas others 
are not so. This result is different from that in [7j. Extending the anlyais to 3cr C.L., all 
cases are consistent with the current fit data. We have presented the numerical results for 
the predictions of 5d in terms of either s\ 2 or for those cases. 
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